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We propose a family of local VT- symmetric photonic lattices with transverse index gradient w, 
where the emergence of stable Bloch-Zener oscillations are controlled by the degree of non-Hermiticity 
7 of the lattice. In the exact VT -symmetric phase we identify a condition between oj and 7 for 
which a wavepacket self -imaging together with a cascade of splittings and giant recombinations 
occurs at various propagation distances. The giant wavepacket recombination is further enhanced 
by introducing local impurities. 

PACS numbers: 


Introduction - Non-Hermitian wave physics and specif¬ 
ically its parity-time (VT) symmetric ramifications [I], 
has attracted a lot of attention in recent years. The 
main observation was that a non-Hermitian Hamiltonian 
PL that commutes with the joint TT-symmetric operator 
may possess an entirely real spectrum. Specifically it was 
shown that below a critical value UvTi °f the parameter 
7 controlling the non-Hermiticity of PL, the spectrum is 
real and the eigenfunctions of PL are eigenfunctions of 
the VT- symmetric operator. In the opposite limit the 
spectrum becomes partially or completely complex while 
the eigenfunctions cease to be eigenfunctions of the VT 
operator. The first domain was coined the exact VT- 
symmetric phase while the latter was coined the broken 
TT-symmetric phase. The transition point 7 = Uvt is 
known as an exceptional point (EP) singularity where 
both the eigenfunctions and eigenvalues coalesce. 

The impact of these ideas is well documented in var¬ 
ious physical settings ranging from matter waves El Ei 
and magnonics [3j to optics [ME], electronics jT9j and 
acoustics [2D]. In fact optics and electronics, where VT- 
symmetric set-ups can be realized by judiciously balanc¬ 
ing gain and loss regions of a system, have provided an ex¬ 
cellent playground for experimentally testing many theo¬ 
retical ideas [MHOS- Among these theoretical predic¬ 
tions jlT, and subsequent experimental realizations 
was a new type of Bloch Oscillations which were unsta¬ 
ble. They either amplified or attenuated since the prop¬ 
agating constants in the associated PT-symmetric lat¬ 
tices became immediately complex (the system entered 
the broken PT-symmetric phase) once a transverse index 
gradient was introduced. 

In this Letter we introduce a class of photonic lattices, 
whose building blocks are PT-symmetric dimers with a 
transverse index gradient ui (see Fig. 0 )- These (quasi- 
one-dimensional) lattices respect a local V,iT -symmetry 
associated with each individual dimer. Despite the lack of 
global PT-symmetry they still have parameter domains 
for which their eigenvalues are real i.e. they are in the ex¬ 
act PT-symmetric phase. In this domain they support 
a new class of stable PT-symmetric Bloch-Zener oscil¬ 
lations which, allow for periodic wavepacket self-imaging 
whenever the choice of the w —7 parameters impose a syn¬ 


chronous behavior between the Zener tunneling and the 
period of Bloch-Oscillations. These synchronous Bloch- 
Zener oscillations experience a cascade of splittings and 
giant beam recombinations which are further enhanced 
in the presence of localized defects. 

Theoretical Model- We consider the photonic lattice 
of Fig. 0 - An experimental implementation of the in¬ 
dex gradient for such a set-up has been realized in Ref. 
[2T1122] . Each waveguide supports only one propagating 
mode, while light is transferred between waveguides via 
evanescent tunneling. The connectivity of the array is 
such that each amplifying (dissipative) waveguide of a 
dimer is coupled, with a coupling constant y, to both of 
the adjacent dimers’ dissipating (amplifying) waveguide. 
In addition we assume an intra-dimer coupling a. We 
will assume that a > A. The diffraction dynamics of the 
evolving electric field amplitude H ! n (z) = {a n {z),b n {z)) T 
of the n th dimer along the propagation direction z, in the 
paraxial description, satisfies the following Schrodinger- 
like equation 

.dan + _ ^ an + Qp) n + 2(6„_i + b n+ 1 ) = 0 

dz 2 /-|\ 

■db n , . „ A . n 1 ’ 

i -1- (nu + iu)b n +aa n + — (a n _i + a n+1 ) = 0 

dz 2 


where a n (b n ) is the field amplitude at the gain (loss) site 
of the dimer. Although the system described by Eq. 0 
does not respect a global TT-symmetry (due to the index 
gradient), nevertheless there is a local VdT symmetry 
that it is satisfied by each individual dimer. 

Spectral Analysis- It is instructive to start by studying 
the dispersion relation of the system in the absence of the 
transverse index gradient i.e. w = 0. Using the Fourier 
transformation a n {z) = f_ n a q (z)e' lqn dq (similarly 
for b n (z)) Eq. ([I]) takes the form: 


,d_ ( a q (z) 
1 dz V b q (z) 


*7 ~v q \ ( d q (z) \ 

-Vq -i't ) \ b q {z) ) 


( 2 ) 


where v q = a + Acos(q). The dispersion relation £ a (q) 
(longitudinal propagation constants) is obtained by cal¬ 
culating the eigenvalues of the 2x2 matrix in Eq. 0: 

£<j{q) = <J\/(a + Acos[q \) 2 - 7 2 (3) 
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FIG. 1: A) The photonic lattice with local VT - symmetry. B) 
The associated dispersion relation for 7 = 0 (blue), 0 < 7 < 
7 PT (red), 7 = 7 P 7 ~(black), and the solid/dashed green lines 
are the real/ imaginary part for 7 > r )' PT . 


where q € (—7r, 7 t] and a = ± indicates the upper/lower 
band. For 7 = 0 the minimal spacing between the two 
bands S = 2 (a — A) occurs at q = ± 7 r. As 7 increases the 
minimal band separation shrinks until the edges touch 
at 7 —> = a — A where an EP degeneracy occurs 

(see Fig. 00 - For 7 > 'y-pT we enter the broken VT- 
symmetric phase and the eigenvalues appear in complex 
conjugate pairs. Below we will focus our analysis on the 
parameter domain for which the spectrum is real (exact 
PT-symmetric phase). In this domain, the eigenvectors 
| rr) of the 2x2 matrix of Eq. ([ 2 ]) take the form 


r) = 


V 2 cos 9 


o —icr0/2 

ae i<?8/2 


0 = arcsin — 


(4) 


and they are also eigenvectors of the VT -operator [Tj. 

When ui 0 the two bands are replaced by two in¬ 
terleaving Wannier-Stark ladders £± = £± + nuj where 
n = 0, ±1, • • •. The offsets £q determine the relative en¬ 
ergy distance between the two ladders and can be eval¬ 
uated numerically from a direct diagonalization of the 
effective Hamiltonian H that describes the paraxial prop¬ 
agation of our system Eq. ( 0 . Moreover, in contrast to 
the u> = 0 case, the system possess multiple exceptional 
points. 

Let us look at the case A = 0. In this case the longitu¬ 
dinal propagation constants £„ are organized in doublets 
associated with the n th isolated dimer: 


+ n = ± sJoA -i 1 (5) 


degeneracies with more remote dimers occur at 7 = 0 
while the previous ones evolve towards larger values of 
7 . The index m, defining the number of degeneracies for 
A = 0, will be used later on in order to delineate the w —A 
parameter space of our system, Eq. (0, into domains of 
broken VT -symmetry (i.e. number of instability regions) 
occurring as 7 increases. 

In Fig. we present a density plot for 7 ™“, asso¬ 
ciated with the first EP, versus uj and A. The purple 
horizontal lines indicate the w m -values discussed previ¬ 
ously. For each such domain, we plot in Figs. 03-05, a 
typical spectral behavior (for fixed A, uj) of the eigenen- 
ergies of the system Eq. (0 versus 7 . We see that the 
number of instability regions is described by the index 
m. 




FIG. 2: A) Numerical results for the first EP of the system 
for a = 1, versus uj and A. B-E) show numerical results for 
the real(black) and imaginary (blue) spectra in domains 1 , 2 , 
3, and 4: for A = 0.5 and a = 1. In B) uj = 2.5, while in C) 
uj = 1.5. In D) and E) w = 0.85 and u> = 0.55 respectively. 
The sub figures in A) are plots of Eq. ( |5|6| ) respectively for 
the same parameters as B) and C). The purple lines have been 
added in A) to illustrate the instability bubbles for small 7 at 
the domain borders 

The domain m = 1, associated with one instability re¬ 
gion (see Fig. 03), can be understood within the frame¬ 
work of a single dimer, see Eq. (0. The latter is also 
plotted at the inset in Fig. 0V. Domain m = 2, can be 
analyzed using two coupled dimers subjected to a gradi¬ 
ent uj: 


For uj > 2a the spectrum is non-degenerate for any value 
of 7 7 ^ a (for 7 = awe have multiple EP degeneracies). 
However for w = 2a we have a degeneracy at 7 = 0, 
where £„ = £^±1 ■ Furthermore for uj = a another (sim¬ 
ple) degeneracy develops at 7 = 0 where now = £„± 2 - 
At the same time the previous degeneracy at 7 = 0 for 
u = 2a, ’’evolves” towards 7 = It is straightfor¬ 

ward to show that for uj m = 2s, where m = 1,2,3,..., 


a 0 f 

nuj + *7 =4 0 

Y (n — 1 )uj — *7 a 

0 a (n — l)w + i'y 

Direct diagonalization of the above Hamiltonian gives: 

£ n = £0 + nljJ \ £q = ±VX ± Y + \jj ( 6 ) 


H = - 


noj — 27 
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where X = (^) 2 + a 2 — 7 2 + (f) 2 , and Y = 

^ A 2 a 2 + (a 2 — 7 2 )w 2 . Eq. (I6| is plotted in the inset 
of Fig. 0 and it describes qualitatively the features (i.e. 
two instability domains) shown in Fig. |2j: associated with 
the system Eq. 0 . Other domains m = 3,4, • • • can 
be explained by analyzing a system of three, four, etc. 
coupled dimers. Below we will concentrate only in the 
parameter space for which the system is in the exact 'P'T- 
symmetric phase (stable domains). We point out again 
that this feature is absent in Ref. m where the system 
is always in the broken VT -symmetric phase. 

Dynamics- To study the dynamics, we have numer¬ 
ically simulated the propagation of a broad Gaussian 
beam for different values of 7 < 'y-pT- We have assumed 
a normal incident, so that at the input plane 2 = 0 the 
beam has excited mainly the first band in a spectral in¬ 
terval around go ~ 0. We first consider the case of 7 = 0 
where the band-gap S = 2 (ct — A) is large enough to allow 
us to neglect Zener tunneling (ZT). According to the ac¬ 
celeration theorem, the transverse propagation constant 
q increases up to ±7r where the wavelength satisfies the 
Bragg condition associated with the underlying periodic 
potential. The wave is then Bragg reflected at propaga¬ 
tion distance 2 = tt/lo and travels in the opposite trans¬ 
verse direction toward lower index sites where it experi¬ 
ences a total internal reflection. The process repeats itself 
leading to a periodic motion which can be considered the 
optical analogue of Bloch Oscillations. The oscillation 
period can be easily estimated using the above consider¬ 
ations and it is zg = 2 tt/uj. The above qualitative picture 
is nicely reproduced in Fig. UK for 7 = 0 and w = 0.231. 

As 7 increases the band-gap <5 becomes smaller and ZT 
between the two bands at their edges q = ± 7 r cannot be 
neglected any more. The associated spreading scenario 
is depicted in Figs. [3)3, C for ui = 0.231 and two different 
values of the gain/loss parameter 7 = 0.405 and 0.443. 
In this case the beam will experience a ZT at distances 
2 ^ = (2n + l)7r/w, where n = 0,1, • ■ ■. Let us discuss 
in more detail the first ZT event at z^ = tt/co. For dis¬ 
tances 2 < z'z ^ the beam is mainly trapped in the lower 
band and propagates along the direction of the local gra¬ 
dient d£~/dq. At 2 IK, due to the tunneling, the beam 
splits into two beams one characterized by the lower band 
and the other by the upper band. While the beam associ¬ 
ated with the lower band reverses direction via Bragg re¬ 
flection, the beam associated with the upper band follows 
a parallel trajectory with d£ + /dq. These two beams will 
again change direction due to total internal and Bragg 
reflections respectively. They recombine at the second 
tunneling point at distance z^ = 37 t/ui. The recombi¬ 
nation process is more complicated as now both occupied 
bands experience coherent interference. We have found 
that at some distances zr = z^ (marked by the third 
green line in Fig. UK) these recombinations can lead to 
a giant power focus (the total power is plotted with red 
line in the z-axis of all upper Figs. [3]). The superposition 
of ZT with Block Oscillations can, in general, result in 


an asynchronous process which destroys exact revivals of 
the initial packet. Nevertheless we find that wavepacket 
self-imaging is achieved for some values of u> — 7. This is 
the case for example for the parameters used in Fig. HP 
(see distance 2 si indicated by orange line) as opposed to 
the results shown in Fig. [3)3 where the self-imaging is 
not observed. 

The dynamics is best analyzed in terms of the Floquet- 
Bloch (FB) eigenvectors of the effective non-Hermitian 
Hamiltonian H that describes the paraxial evolution of 
our system Eq. 0 - Let us indicate, using Dirac’s no¬ 
tation, the FB modes associated with the propagation 
constant (eigenstate of H) £° as \£°). They constitute 
a bi-orthogonal basis and satisfy the following relations 
which are dictated by the symmetric nature of H 




< 7 * 

n 


00 

£m) = £ £ \£°)(£n\ 

(j—zh n= — 00 


1 , ( 7 ) 


where * denotes complex conjugation. Moreover, it is 
easy to show that the FB modes, in the position space 
representation satisfy the following periodicity relation 

{»,l + k\E^ k )=(n,l\£Z) (8) 

where {|/i, ?)} is an orthonormal basis defined by two 
indexes (fJ.,1) with the first index representing the ‘gain’ 
(/r = 1 ) or ‘loss’ (n = 2 ) waveguide while the second one 
denoting the label for the dimer. 

Next, we expand the initial preparation |(0)) 
in the FB basis. The expansion reads | \J/ (0)) = 

E CT =± T,n=-oo c n Kn) wh ere < = (£°*\ 'F(O)). Thus the 
evolving beam is 


1 * (*)> = £ £ 


-i£° 


IC>- 


( 9 ) 


cr—± n =—00 


We now project the evolving beam Eq. 0 to the 
Wannier-Bloch basis |<r, g) = \a) (g> |g) where |g) = 
^ 7 = EE -00 10 ellq s P ans the quasi-momentum space: 


(ct ,g| ^( 2 )) = e zt « z {C (w 2 + g) (a , g| £ 0 


( 10 ) 


+ e-< £ « + Whc+( W 2 + g)(<T ,g|£ 0 +)} 

where we have used the notation (cr, g | = (c | <E) (g | , and 
(a\=(VT\a)) T m- Moreover, the coefficients C a (ojz + 
g) — EpE-oo Cpe~ lp( - Ulz+ A satisfy the periodicity relation 
C a (w2 + g + 27 t) = C a (c jjz + g). 

Equations ( 9|10 ) provide an explanation for the recom¬ 
bination and self-imaging events. They indicate that the 
evolving beam is, in general, not periodic as a function 
of the propagation distance 2 and it is characterized by 
two propagation scales: The first one is the Bloch pe¬ 
riod, zb = 77 , originating from the periodicity of the 
C a functions. The second scale 


z E = 


°0 °0 


is asso¬ 


ciated with the minimal energy spacing in-between the 
two Wannier-Stark ladders and arises from the nontriv¬ 
ial relative phase appearing in front of the second term 
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FIG. 3: Here we show various dynamical evolutions of the lattice where a = 1, w = 0.231, and A = 0.6, for an initial Gaussian 
wavefront of the form a„(0) = b„( 0) = e~ n ' 10 . In the upper graphs the x-axis is the site index where the the amplifying site 
is juxtaposed with the attenuating site on the right ( a n (z ), 6 „(z)) for each index n. The red line in the z-axis shows the total 
power of the lattice while the temperature map color of the plot corresponds to the individual site power (|a„(z)| 2 , \b n (z)\ 2 ). 
The green lines mark the first three zz while the orange lines mark the expected self-imaging time zsi. In the lower graphs the 
gray and pink lines correspond to the normalized relative-power on the right and left of the recombination point, n = 7. The 
black and red dots correspond to the upper and lower band-projections, normalized to one at each instant z. In A) 7 = 0, in 
B) 7 = 0.405, while in C) and D) 7 = 0.443. In D) a defect with strength e = 0.25 is included in the dimer with index n = 7 
and centered around z = ^ with a total length of A z = 8 . The presence of the defect disrupts the expected revival at zsi (see 
the lower subfigure C where zsi is marked with orange line). At the same time it results in a huge power recombination (the 
power pick is at least two times bigger than the one shown in case C) at zr (indicated by the third green line). 


on the rhs of Eq. (|10j). There are u> — 7 values for 
which these two propagation-scales are rationally related 
to one another i.e. ze/zb = N/M. This condition leads 
to a self-imaging of the initial preparation at propaga¬ 
tion distances z$i = Mze = Nzb- For instance, when 
M = 1, N = 5 the initial wavepacket is reconstructed at 
the propagation distance zsi = see the orange line 
in Fig. w . Moreover, a giant power focus (third green 
line) occurs at the recombination event which is between 
two successive self-imaging events. 

A deeper insight of the cascade of recombination events 
can be achieved by evaluating the band contributions of 
the evolving beam. Using Eq. (10), the band contribution 
P a (z) can be calculated as 

P°{z)= \(a,q\^(z))\ 2 dq (11) 

J —TV 


The band projections are plotted on the lower row of Fig. 
[ 3 ] as black (upper) and red (lower) points, where we nor¬ 
malize P~(z) + P + {z) = 1. In the same figures the pink 
and grey lines correspond to the relative power (normal¬ 
ized P R {z) + P L (z) = 1) on the left and right of the 
recombination point, dimer index n = 7, as a function 
of propagation z. This experimental observable strongly 
correlates with the band projections where the pink (left) 
and grey (right) lines correspond to the lower and upper 
bands respectively. The distances z = z^ where they 
demonstrate an oscillatory behavior coincide with the po¬ 


sition where Zener inter-band transitions of power occurs 
according to the semi-classical picture of splittings and 
recombinations discussed earlier (see green lines). 

We have also investigated the effect of a localized de¬ 
fect in the creation of these intense recombination points. 
In general, a defect will devalue the maximum of the to¬ 
tal power; however, a strategically placed defect at one 
of the recombination distances z^ can lead to further 
enhancement of the power peak (see third green line at 
zf in Fig. [ 3 J 1 ) and associated power pick). We interpret 
this phenomenon as resulting from quasi-momentum ran¬ 
domization due to the scattering from the defect prior to 
the Bloch-Zener recombination. This leads to a very fo¬ 
cused recombinations with all power concentrated in a 
very narrow lattice domain. 

In conclusion, we have investigated stable Bloch-Zener 
oscillations in a non-Hermitian lattice with local VT- 
symmetry. We have found that an initial beam expe¬ 
riences a cascade of beam splittings and recombinations 
where the re-concentrated power can exceed the initial 
value due to the non-Hermitian nature of the dynamics. 
At the same time we have found that a judicial selection 
of the index gradient to and the gain/loss parameter 7 
can result in perfect self-imaging of the initial packet at 
distances dictated by these two parameters. This plat¬ 
form can open up new possibilities for the realization of 
reconfigurable beam splitters, interferometers and imag¬ 
ing processing. 
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